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Abstract. The relationship between two different asymptotic techniques, namely, the 
Wigner-Weyl kinetic formalism and the complex geometrical optics method, is addressed 
within the framework of semiclassical theory of wave propagation. More specifically, 
in correspondence to appropriate boundary conditions, the solution of the wave kinetic 
' equation, relevant to the Wigner-Weyl formalism, is obtained in terms of the correspond- 

ing solution of the complex geometrical optics equations. In particular, this implies that 
the two considered techniques yield the same wavefield intensity. Such a result is also 
' discussed on the basis of the analytical solution of the wave kinetic equation specific to 

Q ' Gaussian beams of electromagnetic waves propagating in a "lens-like" medium for which 

i the complex geometrical optics solution is already available. 

> 

1. Introduction 
VO ! 

In the framework of the semiclassical theory of waves [1-5], i.e., short wavelength asymp- 
totics, the uniform (global) description of the wavefield is complicated by the formation of 
caustic singularities [1, 2, 6-8]. Although a complete and deep understanding of the wave- 
field structure near caustic regions is obtained on the basis of catastrophe theory El El 
and the unfolding of the corresponding singularities can be treated by means of symplec- 
■ tic techniques 0121111, the application of such methods to realistic cases, e.g., to waves in 

P" 1 ! magnetically confined plasmas , appears rather difficult. Therefore, with specific regard 

to physical applications, several asymptotic methods have been developed which yield nu- 
merically tractable equations, though being limited concerning the global properties of the 
asymptotic solutions. 

Such asymptotic techniques can be classified into two different families, depending on 
whether the relevant wave equation is described in the phase space, microlocal techniques, 
or directly in the configuration space where the wavefield is defined, quasi-optical methods. 

This work aims to give a detailed comparative analysis of two such techniques, namely, 
the Wigner-Weyl kinetic formalism [11-13] and the complex geometrical optics (CGO) 
method [14-17] which can be considered as benchmarks for microlocal and quasi-optical 
methods, respectively. 

Specifically, in Sec. 2, the Wigner-Weyl formalism and the complex geometrical optics 
method are reviewed and compared. In particular, it is pointed out that, within the Wigner- 
Weyl formalism, physically meaningful solutions should have a specific form, referred to as 
momentum distribution, which is characterized in Sec. 3. 

On the basis of the mathematical properties of momentum distributions, our main result 
is obtained in Sec. 4. In particular, it is shown that, in correspondence to appropriate 
boundary conditions, there exists a specific asymptotic solution of the wave kinetic equation 
relevant to the Wigner-Weyl formalism that can be written in terms of the corresponding 
solution of the complex geometrical optics equations. This allows us to relate the two 
considered methods as well as to determine the specific class of boundary conditions for which 
they are equivalent. In Sec. 5, this general result is illustrated by means of an analytically 
tractable example, i.e., the propagation of a Gaussian beam of electromagnetic waves in an 
isotropic "lens-like" medium. In conclusion, a summary of the main results is given in Sec. 6. 
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2. The Wigner-Weyl Formalism versus the CGO method 

In order to set up the framework, let us discuss the relevant boundary value problem for 
a generic scalar pseudodifferential wave equation, together with the required mathematical 
definitions. 

Specifically, we will consider the case of a scalar (real or complex) wavefield ip(x) propa- 
gating in the iV-dimensional linear space K.^ with x = (x 1 , . . . , x ) a generic set of Cartesian 
coordinates and denote by k — (k\, . . . , fejv) the corresponding coordinates in the dual space 
= K . Time-dependent wavefields are included as one of the coordinates can play the 
role of time, e.g., x N — ct, c being a reference speed, and the corresponding dual coordinate 
is related to frequency, e.g., k^ = —lu/c. To some extent, the results for a scalar wave 
equation are valid also for a multi-component wave equation as the latter can be reduced to 
a set of independent scalar equations far from mode conversion regions [jy. 

Thereafter, the Wigner-Weyl formalism will be formulated entirely in the space l w x 
Qgjjvy ^ R 2Ar j with coordinates (x,k), which is viewed as the trivial cotangent bundle 
|2*T] over the configuration space M. N where the wavefield is defined. No explicit reference 
to the propagation direction is made, differently from the classical derivations As 
for the CGO method, it has been originally developed for solving the second-order partial 
differential equation relevant to the propagation of electromagnetic wave beams in stationary 
spatially nondispcrsive media. Hence, we will need to discuss its application to generic 
pseudodifferential equations. In particular, it is shown that the CGO method yields an 
approximation of the wavefield directly in the configuration space, provided that the wave 
equation satisfies an appropriate condition. 

First, let us define the class of wave equations undergone to solution. A pseudodifferential 
wave equation is an equation of the form 

(Dip)(x) = — [ e ik < x - x ' ) d{x,x',k)ip{x l )d N x'd N k = 0, (1) 

which admits propagating wave solutions. The operator D is a pseudodifferential operator 
[3,18-20] acting on the wavefield as a Fourier integral operator (20j characterized by the 
bilinear phase function k ■ (x — x'). Here, d(x, x' , k) belongs to a particular class of smooth 
functions, referred to as symbols, which, roughly speaking, behave like a polynomial in k for 
\k\ large enough. Specifically, a smooth function a(z,k) with z £ R M and k e (M. )' is a 
symbol of order m £ R if for every multi-indices a, (3 there is a constant G a $ > such that 

\d?d%a(z,k)\ < \kL\) m ~ m , uniformly in (z,k) 6 R M x (M. N )', (2) 

and one writes a S 5 m (M M x (E^)'). In virtue of the symbol estimate J2J) the integral in 
||TJ makes sense for -0 £ S'(R N ), the space of tempered distribution |181 121| . Moreover, the 
scale length L characterizes the variations of symbols with respect to the spatial coordinate 
z and it can be eliminated by the rescaling z — » z/L and k — > kL. It is worth noting that 
any linear differential operator with smooth and bounded coefficients is a pseudodifferential 
operator [TT?] . 

Boundary conditions of Cauchy type are given on an (N — l)-dimensional hypersur- 
face E: for simplicity, one can assume S to be the hyperplane {x : x N = 0} where 
the wavefield ipo(x), x — (x 1 , . . . , x N ~~ 1 ), is assigned together with as many derivatives 
ip n {x) — d n ip / d(x N ) n as appropriate. We are interested in semiclassical solutions for which 
only the covectors k with 

k = \kL\ > 1, 

are significant, in the integral in As a consequence, the wavefield should be a highly 
oscillating function on the (large) scale length L and it should correspond to a specific set 
of highly oscillating boundary conditions of the form 

ifo(x) = A (xy So ^, 9 x 5o(x)L~k»1, (3) 
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the amplitude A being slowly varying, that is, \d x A\ ~ |A|/L. 

2.1. The Wigner-Weyl formalism. In the Wigner-Weyl formalism, the Weyl-symbol 
map a w is applied in order to represent the wave equation in x (M. N )' which is 
naturally endowed with a phase-space structure. The Weyl-symbol map transforms an op- 
erator A : S{R N ) -> S'(R N ), S(R N ) being the space of Schwartz's functions [21], into a 
tempered distribution by acting on the Schwartz kernel A(xi,x 2 ) S S'(R 2Ar ) of the operator 
A according to p3 



A( Xl ,x 2 ) i-» cr w (A)(x,k) = J A{x + \s,x-\s)e~ lk - s d N s. (4) 
It is worth noting that for the pseudodiffcrential operator in JTJ the Schwartz kernel is 
D{x u x 2 ) = 7 ^ w f e ik < x ^d{x 1 ,x 2) k)d N k, 



(2nY 

thus, the image of D under the Weyl-symbol map amounts to the formal series of decreasing 
order symbols 

a w (D)(x, k) = J2 l -^9^d(x + \s,x- \s, k) | s=Q , 

a 

where Taylor expansion has been used and dfd^d <E S ,m_ l Q L Series of that kind admit 
always an asymptotic resummation 11H| to a symbol of order m, 

D(x,k) ~ <T W (D)(x,k) e S m , (5) 

which is referred to as Weyl symbol of D. In © the ~ denotes the asymptotic equivalence 
of symbols Ej • On the other hand, one can consider the correlation operator J^j S 
whose Schwartz kernel is given by the tensor product ip(xi)'tp*(x2), then the Weyl-symbol 
map yields the Wigner function 

W{x,k) = a w (S)(x,k) = J ijj(x + \s)il;*(x-\s)e- lk - s d N s. (6) 

The wave equation can be written in the equivalent form DS = and, on applying 
the Weyl-symbol map, one gets Ej 

a w (DS)(x, k)~J2 {2 ^S\ am M D ( X > k)) (figflgW(a, fc)) = 0. (7) 

In the semiclassical limit one has d^d^D — 0(ft m ~'^') so that, on assuming the same 
ordering for the Wigner function the foregoing equation separates into 

D'(x,k)W(x,k) =0, (8a) 

{W(x, k),D'(x, k)} = 2D"(x, k)W(x, fc), (8b) 

where D' e S m and D" s S" 1 ^ 1 are the real and imaginary parts of the Weyl symbol; in 
particular, one has D"/D' = 0(k _1 ) which is the condition for weak absorption and/or 
instabilities ^Sj- Equation l|8a|) is a constraint to Ea. l|8b|) which, on the other hand, has 
the form of a kinetic equation in the x-k phase space, {•, •} being the corresponding Poisson 
brackets. We will refer to the whole system (jHJl as the wave kinetic equation. 

In general, a solution of the wave kinetic equation is a tempered distribution, however, 
one usually restricts the class of solutions to semiclassical measures. This allows us to make 
sense of the integrals of the form Ej 

( A U = TTTTW I ' A{x,k)W{x,k)d N xd N k, 



which expresses the expectation value of a physical quantity represented by the pseudo- 
differential operator A with Weyl symbol A(x,k). In the following we will assume further 



THE WIGNER-WEYL FORMALISM VERSUS COMPLEX GEOMETRICAL OPTICS 



4 



A(x)= I A(x,k)W(x,k) (9a) 



hK*)| 2 = / T^Tyv W(x,k). (9b) 



regularity with respect to x so that the expectation values of physical quantities can be 
defined locally, that is, 

d N k 
(2tt)' 

makes sense as a smooth function in x G M. N . In particular, the wavefield intensity amounts 
to 

d N k 

Such a restriction of the class of solutions is justified as, in general wave propagation prob- 
lems, one needs a description of the space and time profiles of physical quantities. 

Let us now discuss the appropriate boundary conditions for (JHJ) . On the hyperplane E = 
{x : x N = 0} the wavefield ipo(x) has been assigned and one can compute the corresponding 
Wigner function Wo(x,k) — a w (Sq)(x, k), So being the correlation operator associated 
to ipo and k — (fci, . . . , fcjy-i) the coordinates dual to x; then a solution W(x, k) should 
match Wo(x,k) in some appropriate sense. Specifically, one should impose that the local 
value A(x) of any physical quantity evaluated on E is the same whether it is evaluated 
by W(x,k) or by Wo(x,k). Within this formulation, the Weyl symbol A(x,k) should be 
restricted to M 2JV ~ 2 where Wq is defined. One can note that the suitable embeddings T of 
K>2iV-2 into M. 2N such that T lies over S, i.e., 7rr = E with tt : (x, k) >~ > x the canonical 
projection |2J, are of the form T — {(x,k) : x N = 0,kpj = H(x,k)} with H(x,k) a 
generic smooth function; correspondingly, the restriction of a symbol is readily defined as 
A|r(a;, k) = A(x : 0, k, H(x, fc)) . Then the boundary value conditions read 



d N k 



(2tt)' 



k 



-W(x,k)A(x,k) ^ = J (27r)jV-1 W (x,k)A\ r (x,k), (10a) 

which is equivalent to 



W\ s (x,k,k N ) = 2irW {x 1 k) 5{k N - H{x,k)) . (10b) 

The function H cannot be arbitrary as l|10b|) should satisfy the constraint l|8a() which reads 

D'(x,0,k,H(x,k)) = 0, (11) 

hence, the appropriate functions H are obtained on solving the so-called local dispersion re- 
lation evaluated on E. Since E has been assumed to be noncharacteristic, i.e., dD' /dk^ ^ 
on r, the function H is well defined and smooth at least locally, in view of the implicit 
function theorem. On the other hand, it is not unique since (| 1 II) may have multiple solu- 
tions, each one corresponding to a specific branch of the dispersion relation. In virtue of 
the superposition principle for linear wave equations, the total wavefield is a linear super- 
position of the contributions from each branch of the dispersion relation, the coefficients 
being determined by the Cauchy boundary values of the normal derivatives d n ip/d(x N ) n \s. 
Therefore, one has a specific Cauchy boundary value problem for the Wigner function of 
each branch and the sum over all branches yields the total Wigner function. Since one has 

b,b' b b^b' 

where the indices b and b' run over all branches, and, on noting that the average (•) over short 
scale oscillations cancels out the mixed terms, {ip^t/ib) = for b' ^ b, whereas (\^b\ 2 ) = \ipb\ 2 , 
one gets 

dNk u,f u\ f dNk 
^ W b ( X ,k) = J ^-y 



(WaOH = E / 7^ W*(*> fc ) = / T^n W (x, k), (12) 



that is, the projection of the total Wigner function yields the averaged wavefield intensity 
and, thus, it does not account for, e.g., the formation of short-scale diffraction patterns. 
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2.2. The complex geometrical optics method. Let us now turn to the complex geo- 
metrical optics (CGO) method and, in particular, let us discuss its application to pseudo- 
differential wave equations. This is based on approximating the solution of by a smooth 
wave function of the form 

i/j{x) = u(x) e tS{x) = u{x) e -^e iS ^, (13) 

where, according to the semiclassical limit, \d x S(x)L\ ~ k and d x u(x)\ ~ \u{x)\/L. In ad- 
dition to the standard oscillating exponential e lStyX \ the wave object (|13l) exhibits a novel 
scale length w ~ {dx^x)^ 1 which accounts for intermediate-scales variations of the am- 
plitude profile A(x) = u{x)e~^ x ^ with k S> \kw\ ^> 1. In general, such an intermediate 
scalelength w can be determined by both (strong) absorption and diffraction |17|: how- 
ever, in this paper, it is assumed that the medium is weakly nondissipative [cf. comments 
after equations JHJl] so that only diffraction effects are significant. The total short- and 
intermediate-scale variations of the wavefield are accounted for by the complex eikonal func- 
tion S(x) — S(x) + i(f>(x). 

The relevant equations for the three unknown functions u, (f> and S are determined on 
substituting the ansatz l|13|l into the wave equation JJJ . For the specific case for which D is 
a differential operator this is straightforward. On the other hand, for the general case, one 
should deal with the nonlocal response of the operator [22] ■ With this aim it is convenient 
writing (JIJ in the configuration space in terms of the Schwartz kernel, namely, 

T>{x, x')ij){x')d N x' ~ f T)^(^L,x-x')^(x')d N x' = 0, 



where the exact kernel has been replaced by the 2)W — a w ~ 1 {D) where D is the Weyl 
symbol. Actually, one could make use of other symbol maps j^J, yielding asymptotically 
equivalent results; here the Weyl-symbol maps has been chosen for direct comparison with 
the Wigner-Weyl approach. 

Let us further assume that the kernel T)^ (xi, X2) amounts to a distribution smoothly 
dependent on x\ and with compact support in X2 ■ In virtue of the Paley- Wiener-Schwartz 
theorem [21], this is equivalent to assume that the corresponding Weyl symbol D(x,k) 
extends to an entire function of the complex- valued dual vector k — k + ik" , smoothly 
dependent on x G K . From a physical standpoint, the foregoing assumption implies that 
nonlocal effects have a finite range: the response [Dip) (x) of the operator D depends only 
on the value of the wavefield ip in a compact set. 

Within this condition, one can substitute the complex eikonal ansatz <|13fl into Q and 
expand in Taylor series with respect to s = x — x'. As a result one has [22 

D(x,k(x))u(x)-z g- 

i d 2 D(xMx)) ( , i &S(x) &D{xM*)) , ,_ nl -2, 
2 dx*dki HX > 2dx*dxi dhdkj >> 

where k(x) = d x S(x) + id x (f>(x) and, for the Paley- Wiener-Schwartz theorem, the estimate 
\d^D\ = 0{n m ~\P\) is still valid for complex-extended symbols. To leading orders in k, 
and for a weakly dissipative media, i.e., D" / D' = 0(k _1 ), one gets 

D'(x,d x S(x)) =0, (14a) 



dD'(x,d x S(x)) du(x) 
dki dx l 



dh 



u(x), (14b) 



D'(x,k) and D"(x,k) being the real and imaginary parts of D extended in the complex k 
space and, in general, are complex valued. It is worth noting that the foregoing equations 
can be formally obtained from the standard geometrical optics equations, e.g., in the form 
given by Littlcjohn and Flynn by replacing k(x) with k(x). 
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The CGO equations (|14H have been dealt with both by means of the characteristics 
method in the complex domain |15| and on expanding the equations with respect to e = 
|fc"(x)|/|fc(x)| ~ |fc(x)u;| _1 <C 1 17 . In particular, on referring to the latter approach, in 
the weak-diffraction regime e ~ k -1 / 2 , terms up to order e 2 ~ should be considered 
in the CGO equation for the complex eikonal S(x), which, after separating the real and 
imaginary parts, amounts to 

D'(x,h(x)) - \ kUx)h';(x) 32 ° 3 ^ W) = 0, (15a) 

trw £^ra= . " ' (i5b) 

Equations (|15a|l and l)15b(l constitute a set of coupled first-order partial differential equa- 
tions for S(x) and (f>(x) with k(x) = d x S(x) and k"(x) = d x cf)(x). As for the complex 
amplitude u{x), only the lowest order approximation with respect to e is significant, so that 
the real amplitude |u(x)| is decoupled from the phase arg[u(x)] (not considered hereafter) 
and determined by means of the transport equation 



d 

dx l 



dD'(x,k(x)) 
dk, HX)1 



2D"(x,k(x)) \u(x)\ 2 . (15c) 



This is formally the same equation as the geometrical optics transport equation 0, but 
diffraction effects are accounted for through the wavevector-field k{x) which differs from 
that obtained in the geometrical optics. The approximated form <|15l) of the CGO equations 
is the one used in physical applications. Moreover, in the zero-diffraction regime (w > L), 
one has e ~ K , thus terms up to first order only should be considered, with the result that 
equations Ijl5ajl and l|15b|) are decoupled and the whole set of CGO equations fTK|l reduces 
to the standard geometrical optics equations, (j) being effectively zero. 

Equations (|15a|l and l|15b|) are usually solved by computing the characteristic curves 
[181 12T] of l|15a|l with i|15bjl regarded as a constraint with the result that the characteristics 
curves thus obtained resemble the geometrical optics rays Therefore, the appropriate 

boundary conditions should be enough to determine the initial values of the complex vector 
k\^(x) evaluated on the boundary surface E. 

Such conditions are obtained from the Cauchy data © on writing 

Mx) = A (x)e lSo{x) = u {x)e- Mx) e lSo ^ + 0(e), (16) 

for some functions uq and 4>o such that l^uol ~ \uq\/L and l^^ol ~ 4>o/w. From 116JI 
one readily gets the value of the component of the complex vector k = k + ik" tangent 
to S, namely, k(x) = 8 x Sq(x) and k" (x) = d x <j)o(x). The remaining normal component 
is obtained on imposing that the CGO equations (|15|) are satisfied on S; this yields two 
equation for k^ and fc^, viz., 

D'(x,0,k(x),k N )-^ J2 Aj(x,0,k(x),k N )kUx)k' : l(x)=0, (17a) 

i,j<N 

*w0*0 = -J2 M'(*)*iM>M*),fciv), (17b) 

i<N 

where Aij(x, k) E S m ~ 2 and Xi(x, k) E S° are obtained in terms of the first- and second- 
order k derivatives of D' and evaluated at x N = 0. Equation Ijl7a|l is an 0(e 2 ) perturbation 
of the local dispersion relation l|ll|) . hence, it can be solved by 

k N (x)=H(x,k(x))+0(e 2 ), (18) 

and, in correspondence of l|18l) . Eg. (|17b|) yields k'^(x). As in the Wigner-Weyl formalism, if 
multiple solutions are found, one should write the wavefield as a sum of contributions from 
each branch of the local dispersion relation. 
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From the foregoing discussion, one should note that the CGO method yields the solution 
directly in the configuration space, but one should deal with the set of partial differential 
equations (|15J) . the numerical solution of which can be rather cumbersome. Although the 
characteristics technique can be used for Ea. (|15a|) . the constraint (|15bjl should be solved in 
parallel, thus increasing the computational complexity of the problem. As for the global 
properties of the CGO solution, to our knowledge no general result is still available, though 
numerical solutions 17. show that the CGO solution is regular even near focal points where 
the standard geometrical optics solution exhibits a caustic singularity. 

In contrast, the Wigner-Weyl formalism appears better suited for numerical solutions. 
In particular, the wave kinetic equation can be solved along the corresponding Hamiltonian 
orbits in the phase space so that it is reduced to a set of ordinary differential equations that 
require limited computational efforts and the solution thus obtained has a global validity in 
the phase space since the Hamiltonian orbits do not cross each other. In this respect, the 
constraint (|8a|) does not constitute a limitation as D' is a constant of motion. Moreover, 
there is no limitation on the nonlocal response of pseudodifferential operators to which 
the Wigner-Weyl formalism applies. On the other hand, the solution in the phase space 
should be projected into the configuration space and, thus, an integral with respect to the 
momentum k should be carried out numerically. 

Notwithstanding these differences, the Wigner-Weyl kinetic formalism and the complex 
geometrical optics method share a number of features, e.g., the solution of the local disper- 
sion relation 118|) relevant to the CGO method is obtained, to the lowest significant order 
in e, on evaluating the corresponding solution Hll|) for k = k(x). In the following sections, 
it will be proved that one can project the wave kinetic equation from the phase space into 
the configuration space in such a way that the CGO Ea. (|15(l are recovered. 



As discussed in Sec. 2, the solutions of the wave kinetic equation are usually sought in 
the space <S'(]R 2Ar ) of tempered distributions |2U, or in the space of semiclassical measures 
|3J- The first formulation is the more general, whereas the second follows from the physical 
requirement that expectation values (|9a|l are well defined. Moreover, we have pointed out 
that, for general wave propagation problems, physics requires a stronger condition on the 
Wigner function, namely, the expectation values of physical quantities should be locally 
defined according to l)9b|l. For instance, if one deals with a time-dependent wavefield for 
which x N — ct and fcjv — —u)/c, the integral 



yields the wave action density J(x,t) in the space-time |13j . 

In this section, a mathematical characterization of such novel solutions is given and the 
corresponding differential calculus is put forward. 

First, let us note that for any Schwartz function ip(k) £ 5(11^) and for any tempered 
distribution / G S'(R 2W ) one can define a tempered distribution f v £ S'(M. N ) over the 
configuration space only, given by 



where, in general, angle brackets and the integral are alternative ways to denote the action of 
a distribution on the corresponding test function. The distribution / is smooth with respect 
to x £ M. N if and only if f v amounts to a smooth function f v (x). In this case the map 
ip i ► fip(x) for ip £ ^(M^) defines at every point location x £ R N a tempered distribution 



3. A NOVEL CLASS OF SOLUTIONS TO THE WAVE KINETIC EQUATION 





(19) 



f x e S'{R N ) with 



(/*,¥>) = U( x ) 



for <p £ <S(R ). 



(20) 
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This is a consequence of the completeness of S'(M. N ) along with the identity f v (x) — 

lim(/, xs(' — x)<f) where \s is a compact-supported function that approximates the Dirac's 

s — >o 

5-function for 6^0. 

Let us now consider a symbol A e S~°° = n m S m , that is, A fulfills the symbol estimate 
for every order to. It follows that A x = A(x, •) G S(R N ) and for any / G S'{R 2N ) which 
is smooth with respect to x one can define 



and this is a smooth function on M. N as required in (|9 b|) . The definition i|21|) should be 
extended for any symbol A G S m with arbitrary order. With this aim, the class of physically 
admissible solutions [in the sense of l|9bjl ] is restricted. In particular, it is appropriate 
considering the tempered distributions / G ^'(IR 2 ^) that satisfy the following conditions: 

(1) / is smooth with respect to x, 

(2) the restriction f x amounts to a distribution with compact support, i.e., f x G £'(R N ), 
where £'(M. N ) is continuously embedded in 5'(R ) in the weak topology. 

Such a distribution / will be called momentum distribution since for every x its restriction 
f x represents the distribution of momentum k over x. For short let us write / G Ai for the 
space of momentum distributions. 

Within this formulation for every f € M and for every A G S m Eq. (|21(l is well posed and 
defines a smooth function on M N . Let us note that the foregoing definition of the space M. 
is not the optimal one as functions rapidly decreasing in k are also admissible momentum 
distributions. However, in the semiclassical limit these functions can be ignored and only 
compact-supported distributions are significant. 

Let us now address the derivatives of a momentum distribution / G M. First, the 
derivatives with respect to the momentum k are defined throughout every order. Specifically, 
since 



hence, (d%f) x is compactly supported, and, thus, d%f G M. In terms of the integral notation 
the latter result evaluated for symbols reads 



which is the "integration-by-parts" formula. On the other hand, the derivatives with respect 
to x should be dealt with more carefully. For simplicity, let us consider first-order derivative 



hence (df/dx l ) x is compactly supported and df/dx 1 G Ai. The explicit formula for the 
derivative is obtained on noting that for every symbol A G S m , Af G Ai and 




(21) 





(22) 





(23) 
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The same result would be obtained from the definition \im{df /dx l ,xs{- — x)A x ) through 

5 — >0 

straightforward but longer calculations. Higher-order derivatives can be defined by recur- 
rence, but they are not explicitly needed in the following. 

Searching for solutions of the wave kinetic equations ||Hj in the space of momentum 
distributions leads to the weak formulation 



(27T) 



A'* 



D'(x, k)W(x, k) A(x, k) d N k = 0, (24a) 
{W(x, k),D'(x, k)} - 2D"(x, k)W(x, ft)] A(x, k) d N k = 0, (24b) 



with W G M.. Furthermore, we are interested in semiclassical solutions for which only large- 
enough momenta are significant. On recalling that (W x , 1} = (27r) N \i/j(x)\ 2 , we will search 
for solution of (|24Jl in the form 

W(x, k) = (2ir) N f(x, k - d x S(x)) \^(x)\ 2 (25) 

where S(x) and |V'( ;C )| 2 are smooth functions to be determined; in particular, S(x) defines 
a Lagrangian manifold k — d x S(x) in the x-k phase space. Moreover, / £ M. is normalized, 
i.e., J fd N k — (fx A) = lj an d such that 

K a (x)= [ f{x,k)k a d N ~k = 0{w~ M ), (26) 



for any multi-index a. The integrals in H26J1 are well posed since k a = (k — d x S(x)) a 
are symbols of order |a|, and the corresponding quantities K a (x) express the statistical 
moments of the distribution /, i.e., they expresses how important are the deviations k of 
the momentum from the Lagrangian manifold k — d x S(x). In particular, Kq(x) = 1 in 
view of the normalization condition. In the semiclassical limit it is assumed that the scale 
length w characterizing the range of the momentum deviations is large enough as compared 
to Ic^S^a;)!" 1 = |fc(x)| _1 , namely, \k(x)w\ > 1. 

In Appendix, it is proved that, within the weak formulation l|24l) . the momentum distri- 
bution /, satisfying the foregoing conditions, can be represented by the asymptotic series, 
cf. Eq.gnj, 

f{x,k-d x S(x)) ^J2^T~ Kp(x)d%6(k-d x S(x)), (27) 

1 

controlled by the small parameter e = ^(a:)^!" 1 <C 1. It is worth noting that f(x,k — 
d x S(x)j is thus represented by a distribution which is point supported on the Lagrangian 
manifold k = d x S(x) and completely determined by its statistical moments Kp(x). 
In correspondence of the asymptotic expansion (|27(l , Eq. I|24a|l reduces to 



^ 1 d%D'(x,d x S(x)) K a+I3 (x) = 0, (28) 
a! 



as shown in details in Appendix, cf., in particular, Eq. Formally, Eq. (J2HI> constitutes 
an infinite set of algebraic equations for the statistical moments K a (x) characterizing the 
momentum distribution, where each equation, labelled by f3, is expressed as an asymptotic 
series in e; the function S(x) is determined by imposing that the system (|28|l admits non- 
trivial solutions. Equations l|28() are valid for a general momentum distribution / which 
satisfies (|24all . In particular, on setting K a (x) = for a ^ 0, Ea. (|27Jl reduces to 

f(x, ft - d x S(x)) = (2ir) N 5(k - d x S(x)) , 

which is the geometrical-optics-like solution obtained by Bornatici and Kravtsov [23] an d by 
Sparber, Markowich and Mauser [23 , whereas in (|28|l the only nontrivial equation reduces 
to the geometrical optics eikonal equation [S] for S(x), namely, D'[x, d x S(x)) = 0. 
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4. The CGO-like solution of the wave kinetic equation 

On the basis of the asymptotic expansion (|27Jl for a momentum distribution, one can 
prove the main result of this paper, that is, relating the wave kinetic equation to the CGO 
equations for suitable boundary conditions. 

First let us consider the specific momentum distribution for which 

*,(*) = |°' , , a f-M=2n + l(odd), 

KJ \(-l) n (k"{x)) , for \a\ = 2n (even), V ' 

that is, odd-order moments have been set to zero, whereas even-order moments have been 
related to a single vector field k"{x) — d x (f>(x), with </> an unknown smooth function. Cor- 
respondingly, the momentum distribution 1(2711 takes the form 

f(x,k-d x S(x))~ J2 ~^T~ (k"{x)) P dP k 5(k-d x S{x)) 

|/3|=even 



(-iy 



+00 

E <2n)\ 

n=0 v ' 



d 



-1 2n 



S(k~d x S{x)), (30) 



the second identity being obtained by means of the multinomial formula 



1 0! B r J 



^ J\ ai ~n\ 

with a,i = k'/d/dki (no sum over i). Let us note that the momentum distribution (|3(J|) is 
symmetric with respect to the Lagrangian manifold k — d x S(x), as even-order moments 
only appear; in particular, the second order moment —k'l(x)k"(x) for i = j is negative, so 
that such a distribution cannot be interpreted as a probability measure. 

The momentum distribution ((30(1 should be multiplied by VK 3 ')! 2 to get the whole Wigner 
function 1)25(1 . Let us consider the case for which 

|V(*)| 2 = \u{x)\ 2 e- 2 ^\ (32) 

(j)(x) being defined in (|2l)l) and \u(x)\ 2 is ordered according to \d x u\ ~ \u\/L with L ^> 
w ~ l^^l" 1 . As a consequence, w is the shortest scale length characterizing the wavefield 
intensity \ip(x)\ 2 , hence it can be identified with the scale length w defined after Eo. 1(13(1 . 
namely, w ~ w and e ~ e. 

Then one has the following: 



— (aiH hajv) , (31) 



The Wigner function W(x, k) = (2ir) N f (x, k — d x S(x)^\ip(x)\ 2 , with f given by and 
\%l)(x)\ 2 expressed in the form \Hty . satisfies asymptotically the wave kinetic equation in the 
weak formulation with an 0(e) remainder if and only if 

(i) the smooth functions S(x) and 4>(x) satisfy the complex geometrical optics equations 
f75a[) and JTWj) , 

and 

(ii) the smooth slowly varying function |u(x)| 2 satisfies the transport equation \15c]) . 



First, let us prove the statement (i). In view of the ansatz l(29() . all the equations obtained 
from 1(28(1 with (3 such that \(3\ is an even integer reduce to the same equation which reads 



E 

y. I —even 



(-i) J 



dZD'(x,k(x))(k"(x)) a =J2 



(-1)" 
^ (2n)\ 



k'l{x) 



_d_ 

dk. 



2)1 



D'(x,k{x)) = 0, 



(33a) 
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and, analogously, all the equations obtained from (|28() with j3 such that \/3\ is an odd integer 
reduce to 



E 

|a|— odd 



(-i) J 



+ 00 



d%D>(x,k(x))(k"(x)) a = J2 



(-1) 



n+1 



n=0 



(2n + l)! 



2n+l 



D'(x,k(x)) = 0, 
(33b) 

where the second identity in both equations (|33|l follows on using the multinomial formula 
(|31|l . Equations (|33|l constitute a set of two coupled equations for the real functions S(x) 
and <fi(x), and, to lowest significant orders in e, they are the same as the CGO equations 
(I15a|l and (|15b|l : this completes the proof of (i). 

As for (ii), on account of the differential calculus for momentum distributions put forward 
in Sec. 3, the term connected with the Poisson brackets in Ea. (l24b|l should be written in the 
form, cf. Eq.|13), 



r1 N k r) 



(2tt) n dh 



d N k 



ir 



d fdD' 



A - 



d N k dW dD' 



(2ir) N dx 1 \ dk 

Using the specific momentum distribution (|30|l for which 

W(x, k) = (2Tr) N 6(k - d x S{x)) \ip(x)\ 2 + 0( 

yields 



A. (34) 



r! N k 8 



dD'(x,k(x)) 



\iP(x)\ 2 A(x,k(x)) 



dD'(x,k(x)) dA(x,k(x)) 



(2tt) n dh dx % 
?) (35) 
dD'(x,k(x)) dA(x,k{x)) 



_d_ 

dx 1 



dx l 

dD'(x,k(x)) 
dh 



dh 



dh 

mx)\ 2 + 0(e 2 ) 



dx 1 



|V(*)| 2 



A(x,k{x)) +0(e). 



(36) 



The last identity follows on noting that taking the derivative of l|33a|) with respect to x^ 
yields 

dD'(x,k(x)) dkjjx) _ dD'{x.k{x)) dh(x) _ dD'(x,k(x)) 



dh dx 1 dh dxi 

with kj(x) = dS(x) / 'dx 3 . Equation (|36|l implies that 



dx 3 



0(e) 



{W,D'} = (2tt 



d_ 

dx 1 



dD'(x,k(x)) 



dh 



^(x)\' 



6(k - d x S(x)) + 0(e) 



in the weak sense. Hence, from the wave kinetic equation, to lowest order in e, one gets the 
transport equation 



d 

dx 1 



dD'(x,k(x)) 



dh 



2D"(x,k(x)) \*p(x)\ 



which reduces to the CGO transport Ea. H15c|l for |u(x)| 2 , cf. equation Ij32(l . on noting that, 
to lowest significant order, 



d 

dx 1 



dD'(x,k(x)) 



dh 



\u(x)\ 2 e 



2„-20(x) 



-2<p(x) 



d_ 

dx 1 



dD'(x,k(x)) 



dh 



\u(x)\ 



in view of Eq. I|33bl) . This concludes the proof. 

The foregoing result shows that there exists a specific form of the Wigner function for 
which the wave kinetic equation is reduced to the CGO equations. 
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In order to compare the wavefield intensities predicted by the wave kinetic description 
with that obtained on solving the CGO equations, one should complete the foregoing argu- 
ment by discussing Cauchy boundary conditions. With reference to 111 01) and Qlljl one has 
the following: 

Let W(x, k) be the weak solution of the wave kinetic equation JJ£) corresponding to the Cauchy 
boundary conditions 

W\ s {x,k,k N ) = {2Tr) N 5(k - k(x))5(k N - H(x,k(x)))\u (x)\ 2 e- 2 ' 1 ' ^ + 0(e), (37) 

for some smooth functions Sq(x), 4>q(x) > and \uq(x)\ 2 with k(x) = d x So(x) satisfying 
the CGO ordering defined after $16]) . and let S(x), (f>(x) and \u(x)\ 2 be solution of the CGO 
equations $15]) with Cauchy boundary conditions given by the same function So, <fio and \uq\ 2 . 
Then W(x, k) can be approximated according to 

W(x, k) = (2ir) N 6(k - d x S{x)) \u{x)\ 2 e- 2 ^ x) + 0(e) (38) 

in the weak sense of Sec. 3. 

First, let us note that the Cauchy data l|37[l is a particular case of <|10b[) which corresponds 

to 

W (x,k) = {2ii) N - 1 5{k - d x S {x))\u {x)\ 2 e- 2 ^ + 0(e); (39) 

in particular, the Wigner function corresponding to the complex-eikonal wave object i|16fl 
can be written in the form (|39J) . 

In order to prove the foregoing statement, we will make use of the previous result of this 
section. Specifically, we have proved that the Wigner function 

W(x, k) = {2n) N f(x, k - d x S(x)) Kx)| 2 e^\ 

/ being the momentum distribution l|3U|l . solves asymptotically the wave kinetic equation 
in the weak sense within an 0(e) accuracy. Moreover, 

W(x, k) = (2ir) N 5(k - d x S(x)) \u(x)\ 2 e- 24,(x) + 0(e 2 ), 

in view of l|35|) . As for the boundary condition l|37|l one gets 

W\ s (x,k,k N ) = (2Tr) N 5(k~d x S(x,0))5(k N - d x N S(x,Q))\u{x,Q)\ 2 e- 2rt,{xfl) + 0(e 2 ) 

= (2ir) N 6(k - d x S (x))S(k N - k N (x)) \u a (x)\ 2 e- 2 ^ + 0(e 2 ), 

where kpf(x) = d x S(x,0). According to (|18|l . k^(x) = H(x,k(x)) + 0(e 2 ), and, thus, 
<5(fcjv — kjy(x)) — <5(fcjv — H(x,k(x))^j + 0(e 2 ) in the weak sense, so that W(x,k) matches 
the boundary conditions (|37|l . Since the solution of the wave kinetic equation along with 
the Cauchy boundary condition l|3"T|l is unique and since W is an O(e) solution, it follows 
that 

W{x,k) = W(x,k)+0{e), 

which concludes the proof of (|3"51) . 

This implies that, whenever the solutions of both the wave kinetic equation and the 
CGO equations exist, thus, in particular, the Cauchy boundary conditions are of the form 
(|39|) , the Wigner- Weyl formalism and the complex geometrical optics method are equivalent 
within an 0(e) accuracy. In particular, the wavefield intensity predicted by the Wigner- Weyl 
kinetic formalism is the same as that predicted by the CGO method, namely, 

I" d N k 

IV,(X)|2 = J J2^W W{X > k) = l"^)! 26 " 2 ^ + ( 4 °) 

the second identity following from l|38|l . In the next section, an analytically tractable case 
is considered as an example. Specifically, the solution of the wave kinetic equation relevant 
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to the paraxial propagation of a Gaussian wave beam in a "lens-like" medium is obtained 
and shown to be the same as the corresponding CGO solution. 



5. The KINETIC DESCRIPTION OF DIFFRACTION effects for a "lens-like" MEDIUM 
AND ITS ANALOGY WITH THE QUANTUM HARMONIC OSCILLATOR 

Let us address the case of a monochromatic (e~ iwt ) beam of electromagnetic waves 
propagating in a loss less "lens-like" medium |25| with real refractive index n(r, ui) = n(x) = 

no\l — (x/L)~\ 2 . It is assumed that the wavefield is localized near the axis x — of the 
medium, that is, (x/L) -C 1; moreover, the wave electric field is written in the form E(r, ui) — 
y E(x, z; oj), i.e., it is polarized along the y axis and propagates in the x-z plane. The relevant 
wave equation for the wavefield real amplitude E(x 1 z:co) is thus the Hclmholtz equation. 
The corresponding Weyl symbol is real valued and given by 

ui 2 ui 2 r i 

D{x,k x ,k z ) = -{k 2 x + kl) + —n 2 {x)^~{k 2 x + k 2 z ) + —nlyi~(x/L) 2 \, (41) 

thus, the dispersion relation D — yields two branches, to be referred to as the progressive 
and the regressive waves. As for the Cauchy boundary conditions, let us assume that the 
wavefield is purely Gaussian at z = 0, i.e., E(x,0;u>) — uoexp [ — (x — xo) 2 /wq] , wq being 
the initial width, and the propagation occurs along the z axis, so that one should solve the 
dispersion relation D = for k z . On assuming that each branch of the dispersion relation 
carries half of the wavefield intensity, one can consider the progressive wave only which reads 



fc ~ V \LJ Vfc 



where fco = louq/c is the wavevector at x — and the paraxial approximation (fcz/fco) 2 ~ 
(x/L) 2 ~ (w/L) 2 ~A/L<1 has been exploited as relevant to the weak-diffraction regime 
|25j . It is convenient noting that the dispersion relation corresponding to the second form 
of (|42|l can be written as 



j-(k - k z ) - 
k 



which is formally analogous to the dispersion relation relevant to a quantum harmonic 
oscillator |13| with unit mass and 1/feo — > h, 1/L — > luq, loq being the characteristic frequency 
of the oscillator, z — -> t and fco — kz(> 0) — > cu. In particular, the frequency u> corresponds to 
the shifted wavevector fco — k z along the propagation direction z. The shift occurs because 
of the oscillations of the wavefield along the propagation direction z. 

This analogy allows to make use of the well-known solution of the wave kinetic equation 
for the quantum harmonic oscillator |13| to describe the paraxial propagation of a Gaussian 
beam in the "lens-like" medium. More specifically, the solution of the wave kinetic equation 
for the harmonic oscillator corresponding to an initially Gaussian wave packet ip(x,0) = 
(w ^/^/2)- 1/2 exp [ - (x - x ) 2 /w 2 ] is [El 



w 2 (t) = [cos 2 (w i) +e 2 sm 2 {u t)] vJ 2 , (44b) 

where w(t) is the width of the wave packet as a function of time and £ = 2h/muJow 2 n m 
being the mass of the oscillator and xq the initial displacement of the Gaussian from the 
centre of the elastic force acting on the oscillator. 
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Figure 1. The wavefield intensity Ij45(l in the (x, z)-plane as obtained from 
the wave kinetic equations JHJ) (bright regions correspond to high intensity), 
for the case L/zr = 0.5 [cf. equation (j45b|) ] with x — and xq — ^w , 
respectively. One should note that the wave beam exhibits a finite width 
even near focal points (characterized by bright spots) where the geometrical 
optics solution exhibits caustic singularities. 



Correspondingly, the solution for the wave electric field intensity in the "lens-like" 
medium, with the considered launching conditions, is 



(\E(x,z\w)\ J )=i4 



2 W 



w(z) 



exp 



w 2 (z) 



cos 2 (z/L) 



2L \ 2 



sin 2 (z/L) 



^ (x — Xq COs(z/L)Y 

w(z) 2 
1 + ((-Zr' 



- 1 sin 2 (z/L) 



(45a) 



w 2 , (45b) 



with zr = the Rayleigh range in the medium. In Ea. l45b|) . it has been explicitly 

indicated that the solution obtained from the wave kinetic equation amounts to the averaged 
intensity (|E(x, z; oj)\ 2 ), rather than to the exact value |E(x, z; cu)\ 2 since two branches of the 
dispersion relation exist each one carrying half of the wavefield intensity, cf. comments after 
Eas. l|12|l . The intensity (|45a|l and the beam width 145b|) are the same as the corresponding 
quantities obtained from the CGO solution [251 - ^ s a consequence the intensity profile l|45[) 
accounts for diffraction effects as shown in Fig.l. 

One can conclude that, according to results of Sec. 4, the kinetic formalism can be used to 
describe the effects of diffraction on the propagation of wave beams, and, for the case under 
consideration, it yields the same result as the CGO method. Nevertheless, the detailed 
structure of the wavefield, i.e., the oscillations along the propagation direction z, the effects 
of the curvature of phase fronts and the Gouy shift, which are available from the CGO 
solution |25j . cannot be resolved by means of the wave kinetic equation, which instead gives 
the averaged intensity distribution. 



6. Conclusions 

Within the framework of semiclassical wave propagation, two specific asymptotic tech- 
niques have been considered, namely, the Wigner-Weyl kinetic formalism and the complex 
geometrical optics (CGO) method. A detailed comparative analysis of these techniques has 
been given in Sec. 2, for the case of scalar pseudodifferential wave equations, with Cauchy 
boundary conditions. 

In particular, in the Wigner-Weyl formalism, the wavefield is represented in the phase 
space by the Wigner function which is a solution to the wave kinetic equation. In the most 
general case the Wigner function amounts to a tempered distribution. However, physical 
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considerations lead to the definition of a novel class of weak solutions which have been char- 
acterized in Sec. 3. Such specific weak solutions are referred to as momentum distributions 
since, for every point location x in the configuration space, they give the distribution of 
momentum k over x, in the x-k phase space. 

On the other hand, the CGO method yields an asymptotics solution of a pseudodiffer- 
ential wave equation directly in the configuration space, in terms of three smoothly varying 
functions representing, the phase, the wavefield envelope and the amplitude, respectively. 

In Sec. 4, on the basis of the mathematical framework developed in Sec. 3, we have proved 
that, whenever both the solutions of the wave kinetic equation and of the CGO equations 
exist, thus, in particular, the Cauchy data are of the form l|39|) . the former can be ap- 
proximated by a momentum distribution, cf. , equation l|38|) , written in terms of the three 
smoothly varying functions that solve the CGO equations l|15|) . As a consequence, the two 
considered techniques are asymptotically equivalent and, in particular, to lowest significant 
order, the wavefield intensity predicted by the Wigner-Weyl formalism is the same as that 
predicted by the CGO method, cf. Ea. (|4T)|) . 

In addition, one can conclude that the Wigner-Weyl kinetic formalism properly describes 
the wavefield near focal points. This is also shown by comparing the solution of the kinetic 
equation to that of the CGO equations for a specific case, namely, the propagation of 
electromagnetic Gaussian wave beams in an isotropic "lens-like" medium, cf. Sec. 5. In 
particular, the relevant solution of the wave kinetic equation has been obtained on the basis 
of the analogy between the "lens-like" medium and the quantum harmonic oscillator. 
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Appendix A. The asymptotic series expansion of the momentum distribution 

In this appendix the asymptotic series expansion l|27(l of the momentum distribution is 
proved and the corresponding equations l|28(l are derived from the weak form Ij24a(l of the 
dispersion relationship l|8a|l . 

Since both D'(x,k) and A{x,k) are symbols, they are, in particular, smooth functions, 
and one can apply the Taylor's formula 



D' 
A 



(x,k(x))k a 



a = (ai, . 



(x,k(x) + k) = ]T ^ 

| a | < n — 1 

a at) being an iV-dimensional multi-index and 



d%D' 



(x,k) 



(46) 



(x,k) 



-r / dt(l-ty 
a! Jo 



d%D' 
d t A . 



(x, (1 -t)k(x) +tk) 



the remainder of order n relevant to the expansions of D' and A, respectively. More specif- 
ically, on making use of l|46(l to evaluate the left-hand side of l)24a|l . one gets 

f(x, k)D'(x, k(x) + k)A(x, k(x) + k) d N k 

~ E d * D '( x > k< ^)HA{x, k{x)) K a+fj (x) (47) 



a, (3 



where K a (x) = 0(w ) are the statistical moments of the momentum distribution f(x, k), 
cf. Ea. (|26H . In virtue of the symbol estimate (0), symbols are such that, e.g., \d£A(x, k)\ — 
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0(|£;| m ~! a !) in the semiclassical limit \k\ — ► +00 uniformly in x, hence, the asymptotic series 
expansion (|47(l is controlled by the (small) parameter e = \k(x)w\~ 1 . Moreover, on noting 
that 

d%A(x,k{x)) = J d%8(k-k(x)) A(x,k) d N k, 

Eq. (|T7|> takes the form 

D'(x, k)f(x, k - k(x))A(x, k) d N k 



£ d kD'(x,k(x))K a+p (x)d P k S(k - k(x)) 

a, {3 



A(x,k) d N k 



and, in view of the arbitrariness of A[x, k), one gets 

D'(x,k)f(x,k-k(x)) ^E^W- d kD'(x,k(x))K a+0 (x)d^S(k-k(x)), (48) 

in the weak sense. It is worth noting that the derivation of l|48|) does not depend on the 
explicit form of the symbol D'(x, k), thus, on setting D'(x, k) = 1, Ea. (|48|l reduces to 

f(x, k - k(x)) ~ E Kp{x)d?8(k - k{x)) (49) 

which is just the general asymptotic expansion l|27(l of the momentum distribution. 

Going back to Eq. I|24a|l , its solution is obtained on setting the expansion l|48|l to zero and 
exploiting the linear independence of the derivatives of the Dirac's <5-function, thus yielding 
a set of equations for the statistical moments, namely, 

E^ d%D'(x,k(x))K a+f3 (x) = 0, (50) 

a 

which is just Eq. l|28|l . It is worth noting that Ea. (|5U|l can be also obtained on substituting 
(|49|l into H24a|) and exploiting the Leibniz's formula 

d "( D ' A )= E ^L^ D ' d l A 

a+7=/3 1 

which expresses the derivative of a product to any orders. 
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